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1 Introduction 



The spin-statistics theorem has a fundamental role in connecting the quantum mechan- 
ics of many-body systems and quantum statistical mechanics [TJ [2] . The wave function 
of a system of identical integer-spin particles (bosons) is symmetric under exchange, 
while it is antisymmetric for a system of identical half-integer particles (fermions). The 
symmetrization and antisymmetrization of the many-body wave function is directly re- 
lated to the commutation and anticommutation relation of creation and annihilation 
operators in the language of second quantization. Also, the number of particles occupy- 
ing the same quantum state and the statistics are mutually connected. However, there 
is some generalization for the ordinary statistics which is based on the generalization of 
symmetrization of the many-body system wave function [31 0], the Pauli exclusion prin- 
ciple [3 [6], the algebra of commutation and anticommutation relation (TJ El 13 [TUl [TT] , 
and the number of particles occupying the same quantum state [121 13] • There is yet, 
another kind which is based on the generalization of the Boltzmann-Gibbs entropy [H]. 

The thermodynamics of different models can be considered by a qualitative tool, 
namely the thermodynamic curvature, which has been introduced by the theory of 
thermodynamic geometry ^Hl UHl EH HE1 HSl 1201 12H 1221 12HI I2H I2SI 1251 EH I2S]- Ther- 
modynamic geometry can be used as a measure of statistical interaction, e.g. the 
thermodynamic curvature of an ideal classical gas is zero and it is positive and negative 
for boson and fermion gases, respectively, indicating that the statistical interaction of 
these models are non interacting, attractive and repulsive [201 IZU [22]. It has been ar- 
gued that the scalar curvature could be used to show that fermion gases are more stable 
than boson gases and that, therefore, one can utilize the scalar curvature as a stability 
criterion. We notice that the association of the scalar curvature with the stability of the 
system is somewhat tentative and conjectural and proposed by Janyszek and Mrugala 
[2T| [22l [25]. Also, it has been shown that the singular point of the thermodynamic cur- 
vature coincides with the critical phase transition point of some thermodynamic systems 
[T9~l |22"1 [2H [27] . Recently, we worked out the thermodynamic curvature of a thermo- 
dynamic system in which particles obey fractional exclusion statistics [291 EQl EI] . The 
study yielded some interesting results about the thermodynamic properties of these 
systems. In the present paper, we mainly investigate the thermodynamic geometry 
of another kind of intermediate statistic, i.e., g-deformed bosons and fermions. One 
interesting realization of the intermediate statistic is related to the deformation of com- 
mutation and anti-commutation relations, which has led to the g-deformed algebra of 
creation and annihilation operators and is related to the general theory of quantum 
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groups [U [9] . Many authors have considered the statistics and thermodynamics of such 
models [32l [33l EH [35l ESI EZl EHl EHl SQl SH S2] - We construct the thermodynamic 
geometry of g-deformed boson and fermion gases by using the g-calculus introduced by 
Jackson [13] . 

The rest of paper is organized as follows. In Sec. II, we summarize the g-calculus 
and deformed algebra. We derive the internal energy and particle number of g-deformed 
bosons and fermions based on the Jackson derivative (JD) of the logarithm of partition 
function. In Sec. Ill, we construct the thermodynamic geometry of g-deformed bosons 
and fermions. Finally, we consider the singular point of thermodynamic geometry to 
show that the singular point coincides with the critical point, where a phase transition 
such as Bose-Einstein condensation occurs. 

2 Deformed Algebra and g-calculus 

The symmetric g-oscillators algebra is defined by the following equations [32J [33J [31] , 

cd — kqc^c = q~ N (1) 



[c, c] k = [c\c% = 0, 

[N,J] = J, [N,c] = -c, (2) 
Jc=[N], cc f = [1+kN], 

where, [A, 5]^ = AB — xBA, and k — 1 stands for symmetric g-bosons while k = — 1 
for symmetric g-fermions. Also, c, and N are the creation, annihilation, and number 
operator, respectively. In the generalized g-deformation, the g-basic number and the 
Jackson derivative of real functions are defined as 

w - (3) 
^ m= m^. (4) 

Now, consider the following Hamiltonian of non-interacting g-deformed bosons and 
fermions 
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where, is the kinetic energy in each state, fx is the chemical potential and N{ is the 
number operator, which is deformed according to Eqs. (j2j). Mean occupation number 
is defined by 

Applying the cyclic property of the trace and using the g-deformed algebra, one can 
obtain 

[m] = e-*[km + 1], (7) 

where, 77 = pfe—fi). Using algebraic calculations, one can evaluate the mean occupation 
number as [SJ El ES] 

1 / z~ 1 e l3ei — kq~ ks 

n * = =1 ln -1 Re- T~k 

q — q \ z e p 1 — kq K 

Now, we can start with the logarithm of the grand partition function of the g-deformed 
boson and fermion, which is given by 

InZ = -fc^ln(l -kze-to). (9) 

i 

We note that the Hamiltonian is designated to be a linear function of the number 
operator. It is straightforward to show that the thermodynamic quantities in ordinary 
statistics such as boson or fermion gas are given by 

d d 

N = -(^-\nZ) l3 = z(^-\nZ) yi = J2^ (10) 

i 

where, z = e M//fcT = e~ 7 is the fugacity and %ji = e -/3e \ It has been shown that by 
replacing the ordinary derivative with JD, one can find a simple form for the number of 
particles in the g-deformed case [Ml HIJ H2] . In other words, we define the JD derivatives 
with respect to the thermodynamic parameters (/3,j) by 

Vf = -e^f, 2?W = -zVf, (11) 

and immediately the well defined particle number and internal energy is obtained by 

U = -V®\nZ=-k^e i V$]n(l-kze-te)\=Y i eini, 

i i 

N = -V®hi Z = zV^\nZ\ y = ^n h (12) 
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where, n.j has been expressed in Eq. (jSJ). In the thermodynamic limit and in a D 
dimensional box of volume L D with particles satisfying the dispersion relation e = ap a 
where a takes constant values for example o = 2 for non relativistic and a = 1 for ultra 
relativistic particles,the summation can be replaced by the integral. 

. AD roo 

, 13) 

where, A = {2/a) l ^ D ^j^ h is a constant which for simplicity, we will set it practically 
equal to one (A = 1) later. For the symmetric g-deformed statistics, we can work out 
the thermodynamic quantities using the mean occupation number in Eq. (jSJ), 

A D 1 f d (z- x eP*-kq- k \ 

- r(f) q - q -^ H %^ z ^ q) 

A D 1 f°° d_ u _ fz- l e^-kq~ k 



N = i — TiTT / In — — f-r)de 



A r £ g-r 1 ^*'*'^ (14) 



where, 



H n (az, k, q) = Li n (azkq k ) — Li n (azkq k ), (15) 
where, Li n (x) denotes the polylogarithm function. 



3 Thermodynamic curvature of a deformed statistic 
gas 

The geometrical study of the phase space of statistical thermodynamics dates back 
to J. W. Gibbs [15] . The geometrical thermodynamics was developed by Wein- 
hold and Ruppeiner [161 02] ■ Ruppeiner geometry is based on the entropy repre- 
sentation, where we denote the extended set of n + 1 extensive variables of the sys- 
tem by X = (U, N 1 , V, N r ), while Weinhold worked in the energy representa- 
tion in which the extended set of n + 1 extensive variables of system are denoted by 
Y = (S, N 1 , ...,V, N r ) [18J. It should be noted that we can work in any thermo- 
dynamic potential representation that is the Legendre transform of the entropy or the 
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internal energy. The metric of this representation may be the second derivative of the 
thermodynamic potential with respect to intensive variables; for example, the thermo- 
dynamic potential <3> which is defined as, 

$ = $({F*}) (16) 

where, F = (1/T, —^/T, ...,P/T, — ///T). $ is the Legendre transform of entropy 
with respect to the extensive parameter X\ 

' 9S (IT) 



OX 1 ' 

The metric in this representation is given by 

<9 2 $ 

Gl > = dFWFi' (18) 
Janyszek and Mrugala used the partition function to introduce the metric geometry of 
the parameter space [2T] . 

G- = ^ (19) 

where, /3* = F % jk and Z is the partition function. 

According to Eqs. (JHJ), the parameter space of an ideal g-deformed gas is (/?, 7), 
where (3 = 1/kT and 7 = —fi/kT . For computing the thermodynamic metric, we 
select one of the extended variables as the constant system scale. We will implicitly 
pick volume in working with the grand canonical ensemble [21J. We can evaluate the 
metric elements by the metric definition in Eq. ( fl9|) . But, there is a subtle point to 
notice. According to the definition of thermodynamic quantities such as the particle 
number and internal energy, we have to replace the ordinary derivative by JD to obtain 
the correct g-deformed quantities. Therefore, naturally we define the metric element 
by the second JD of the logarithm of the partition function. Also, the first JD of the 
logarithm of the partition function yields internal energy and particle number. Thus, 
the metric elements of the thermodynamic space of an ideal g-deformed gas are given 
by 

G„ = vf{vfh,Z) = -vfu = y £^n i 

i 

G M = vf{VflnZ) = -VfN = Y,^S (20) 

i 

G 77 = V^\V^\xiZ) = -VfN = z^V®m. 
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Replacing the summation by the integral according to Eq. ( !T3|) and evaluating the JD 
of mean occupation number in Eq. (jSJ), we will have the following relations: 



r(^ + 2) ( \ 

G W = F %) ( g _ g -l)2 +3^' fc ' g) ~ g f +3^g >MJ 

r(^ + i) /j-cf+i) / _ 

G /37 = ^ _ g -i) 2 [H^ +2 (zq, k, q) - H^ +2 (zq ,k,q) 

r(£) / _ \ 



(21) 



We have evaluated the metric elements of the two-dimensional thermodynamic ge- 
ometry of a g-deformed gas. For evaluating the Ricci scalar (R) of this geometry, one 
can use the well-known definition of the Riemann (RL k ) and Ricci (Rij) tensors, 

pi a pi a pi _i_ rmrl _ pm-pi 

^ijk ~ u j L ki u k L ji i 1 hi 1 jm 1 ji L kmi 

R ij = R ikj> (22) 

R = G Rij , 

where, is the Christoffel symbol defined by 

^ijk = ^ V^ij,k "I" GikJ G jk,i) j (23) 

where, = dkGij, and the two-dimensional parameter space is defined by (3 and 7. 
It should be noted that we do not use Jackson derivative for constructing a geometrical 
picture; however, we constructed the metric elements by Jackson derivative of the par- 
tition function. The sign convention for R is arbitrary; so R may be either negative or 
positive for any case. Our selected sign convention is the same as that of Janyszek and 
Mrugala, but different from [18]. We have used MAPLE for calculating the thermody- 
namic curvature. The results are provided in the relevant diagrams. Figure [1] shows 
the thermodynamic curvature of the symmetric three-dimensional non-relativistic ideal 
g-deformed boson gas (k — 1) as a function of fugacity in isotherm processes for some 
value of the deformation parameter. It is observed that the thermodynamic curvature is 
positive in the full physical range for an arbitrary value of deformation parameter. An 
interesting phenomenon is observed in g-deformed bosons at a specified value of fugacity 
so that we can argue that there is an upper bound on the fugacity of g-deformed boson 
gas where the thermodynamic curvature is singular. Also, it is obvious from Fig. [2] that 
the thermodynamic curvature has the symmetric property (g -H- g _1 ) and that, for a 
specified value of fugacity, we can find two different deformed boson gases with q < 1 
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Figure 1: The thermodynamic curvature of a 3D symmetric g-deformed boson gas with 
respect to fugacity for an isotherm (/3 = 1) and some values deformation parameter (q = 
0.5 = 2)(blue solid line), (q = 0.7 = 1/0.7) (red dash-dot line) and (q = y/o7? = 
(green dot line). 

and q > 1 but with the same value of thermodynamic curvature. We note that the value 
of thermodynamic curvature for any deformation parameter increases with increasing 
fugacity and that specially for the boson gas (q — 1), this value goes to infinity at the 
related upper bound of fugacity (z — 1). 

Based on the interpretation of the thermodynamic curvature as a measure of sta- 
tistical interaction, we can argue the g-deformed boson gas, as in the case of the boson 
gas, has an attractive statistical interaction. The stability interpretation of the thermo- 
dynamic curvature which is only a conjecture, predicts that the g-deformed boson gas 
may be more stable for larger (smaller) values of the deformation parameter for q < 1 
(g>l). 

Another important point concerns the singular point of thermodynamic curvature. 
We can deduce that for any value of the deformation parameter there is a critical 



8 




Figure 2: The thermodynamic curvature of a 3D symmetric g-deformed boson gas with 
respect to deformation parameter for an isotherm (/3 = 1) and some values of fugacity 
(z = 0.05) (red solid line), (z = 0.1) (green dotted line), (z = 0.3) (blue dash-double 
dotted line) and z = 0.5 (purple dash-dotted line). It is obvious that the related 
deformation parameter at singular point for any value of fugacity is given by q = ^Jz 
for (q < 1) and q = l/y/z where (q > 1). 

fugacity z q where, R is singular. This critical value of fugacity is given by 

* = ( i q < \ < 24 ) 

[ q q > 1 

Similar calculations have been performed for the standard g-deformed boson gas, 
which is defined in jlQ] and is known in the literature, yielding results quite similar to 
those for the symmetric one. The only difference between these two g-deformed bosons 
is in their critical fugacity. In the standard case, the following result can be obtained 
from applying our proposed method or others [3D] . 

z q =\ _ 2 * , (25) 
\ q q>l 
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Figure 3: The thermodynamic curvature of a 3D symmetric g-deformed fermion gas 
with respect to fugacity for an isotherm (/3 = 1) and some values deformation parameter 
(q = 0.3 = l/0.3)(blue solid line), (q = 0.5 = 2) (red dash-dot line) and (q = 0.8 = 1.25) 
(green dot line). 

By setting (k = —1), we can investigate the thermodynamic curvature of g-deformed 
fermion gas. Again, we restrict ourself to the three-dimensional non-relativist ideal gas. 
The results are given in Fig. [3J The thermodynamic curvature of a g-deformed fermion 
gas is negative for any arbitrary value of deformation parameter in the full physical 
range, which indicates that the statistical interaction of the g-deformed fermion gas is 
repulsive as is the case with the fermion gas. Also, based on the stability interpretation, 
we conclude that the g-deformed fermion gas may be more stable than the g-deformed 
boson for an arbitrary value of g. We have calculated the thermodynamic geometry for 
g-deformed boson and fermion gases in other dimensions and in the ultra-relativistic 
limit. Almost similar results have been obtained (not shown here to save space). 

We started with Jackson derivatives (JD) to construct the metric elements. If one 
insists on replacing JD with the ordinary derivative, a new q-deformed geometry may 
then be developed by using JD to obtain a q-deformed Ricci scalar. Jackson derivative 
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of the metric elements are given in the following equations for an ideal g-deformed boson 

gas, 

_|_ 3) /j-(^+ 3 ) / \ 

Gppj = °, D . 7 zr [IHd +4 (z,k,q) - H R 4 (zq 2 ,k,q) - Hd A {zq' 2 ,k,q)) , 

1 ( -j) \Q ~ 1 J v " CT CT / 

G /3/?,7 = p/gx ^ _ ( 2 #f + 3 (^ k > O) ~ HR +3 (zq 2 , k, q) - H^ +3 (zq' 2 , k, q)J , 

G 77,/3 = ifzm 7 ( 2 HD 2 (z,k,q) - HD 2 (zq 2 ,k,q) - H R2 (zq 2 ,k,q)\ , 

1 ("2 j W ~~ 9 J V CT CT 7 

G 77,7 = wd| / ~ — =T\3 (^^o+ifoM) ~ H & +l (zq 2 ,k,q) - H o +l (zq' 2 ,k,q)) . 

f ("2 ) I? — 1 ) V CT CT J 

(26) 

It is straightforward to calculate a q-deformed thermodynamic curvature. This q- 
deformed Ricci scalar has properties similar to the usual scalar curvature such as posi- 
tivity for g-deformed bosons and negativity for g-deformed fermions. The only difference 
will be in the upper bound of fugacity, which will be 



■ q > 1 



Although a g-deformed geometry seems to be interesting from theoretical and mathe- 
matical viewpoints, no appropriate physical model is available. 



4 Singular point and condensation 

It has been shown that the singular point of thermodynamic geometry coincides with 
the phase transition in several thermodynamic systems. For example, the thermody- 
namic geometry of Van der Walls gas and that of the Ising model are singular at phase 
transition points. Also, for an ideal boson gas, the Bose-Einstein condensation occurs 
at z = 1, where the thermodynamic curvature is singular. Recently, we have shown 
that for a kind of fractional exclusion statistics, the singular point of thermodynamic 
curvature corresponds to the condensation of a non-pure bosonic system [31] . Now, we 
investigate the singular point of the g-boson gas. It has already been shown HI] that 
the definition of heat capacity leads to an upper bound on the fugacity of g-deformed 
boson gas, while there is no bound on the fugacity of g-deformed fermion gas. We expect 
a condensation to appear where the upper bound of fugacity for a g-deformed boson 
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gas coincides with the singular point of the thermodynamic curvature (Eq. ff24|) ). In 
this formulation, it is straightforward to obtain the phase transition temperature. We 
restrict ourselves to the q < 1 cases, which is simply extendable to the q > 1. According 
to the particle number in Eq. ( TT^l) . the critical temperature T c at the phase transition 
point is given by, 

-| cr/D 

(28) 



k B T.= &r' D "' 



7T 



tr/2 



nT(D/2) q-q 



-i 



T(D/a) Ho +1 (q*,l,q) 



where, n = N/L D is the density of particles and assumed to be constant. It has been 
shown that there is a condensation for the boson gas at a finite temperature in the three 
dimensional space with non-relativistic particles, while the phase transition point goes 
to zero in the two-dimensional space. Generally, Bose-Einstein condensation does not 
occur at D/a < 1 for an ideal boson gas. We investigate the possibility of condensation 
for an ideal g-deformed boson gas in various cases. The important parts of Eq. ( 128]) . 
which shows the dependency of temperature on dimension, are the dispersion relation 
and the deformation parameter: 

A ^ a ' 9) = / 7n S = T- FwK-- Tv ( 29 ) 

HD +l (q 2 ,l,q) LiD +l (q 3 ) - LtD +1 (q) 

<r a a 

It is interesting that where q ^ 1, A(D,a,q) has a finite value for any dimension. In 
the boson gas limit, this function reduces to 

)hnA{D,a,q) = —± m , (30) 

a 

which is zero for D/a < 1. In other words, for general values of the deformation 
parameter, one can find a condensation in all dimensions with non-relativistic or ultra- 
relativistic particles, while for the boson gas (q — 1), the condensation occurs in three 
(and higher) dimensions and for ultra-relativistic particles, in two dimensions [HI B5j 
H6]. Figure (j3J) demonstrates the critical temperature as a function of deformation 
parameter in various dimensions for non-relativistic particles. One can find that, in all 
dimensions, there exists a finite critical temperature for a q-deformed boson gas. The 
exception is one and two dimensional cases in the limit of boson gas (q = 1). Also, it 
is obvious that the finite critical temperature reduces when the deformation parameter 
increases and that in some cases, it goes to zero, which excludes condensation in such 
cases. 

According to the general form of the metric elements and their derivatives, the 
thermodynamic curvature of deformed boson and fermion gases is a function of and 
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Figure 4: Critical temperature as function of deformation parameter in D — 1 (blue 
dashed line), D = 2 (purple solid line), D = 3 (dash-dotted line) and D = 4 (dotted 
line) for non relativistic particles. B = ^ 2 /g^ fe2 is a constant. 

z in the following general form 

R = /3%f(z,q,D,a). (31) 

Therefor, in the low temperature limit and at a fixed fugacity, the thermodynamic 
curvature diverge as R ~ T~%. This is a general behavior for quantum gases such as 
boson, fermion, fractional exclusion statistics, Gentile statistics, as well as g-deformed 
boson and fermion gases. The sign of the thermodynamic curvature in this limit depends 
on the statistical interaction. In our sign convention, it is positive for boson and q- 
deformed boson gas while it is negative for the others. 

5 Conclusion 

The study of thermodynamic geometry of thermodynamical systems has an archaism as 
long as three decade. The investigation of known thermodynamic systems has created 
some application for a new quantity, namely thermodynamic curvature. The thermo- 
dynamic curvature is served as a measure of statistical interaction, a sign for phase 
transition, a tentative measure of stability and in general the critical properties. In this 
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paper, we have endeavored to consider the thermodynamic geometry of symmetric q- 
deformed bosons and fermions. Therefore, we can argue that the statistical interaction 
of g-deformed boson gas is attractive, while it is repulsive for g-deformed fermion gas. 
Also, due to the stability interpretation of thermodynamic curvature, the g-deformed 
boson gas for large value of deformation parameter (q < 1) might be more stable than 
the small value. It means that an ideal boson gas is more stable than the deformed 
one. Specifically, up to the validity of the stability interpretation, we argue that the 
g-deformed fermion gas could be more stable than the g-deformed boson gas for any 
arbitrary value of deformation parameter. Also, investigation of singular point of ther- 
modynamic curvature represent that the critical fugacity for any arbitrary value of 
deformation is coincide to the upper bound of fugacity, which is given by the heat 
capacity definition. Therefore, we can define the phase transition temperature such 
as Bose-Einstein condensation temperature. We showed that g-deformed boson gas 
(g 7^ 1), may condensate in finite temperature in all dimension with particles satisfying 
the arbitrary dispersion relation. The existence of the finite phase transition tempera- 
ture for g-deformed boson gas even in low dimension is an interesting and new result of 
present paper. 

Obtaining the critical temperature for standard g-deformed boson gas [ID] is straight- 
forward such as previous one and we elude to calculate and represent it because of entire 
similarity. 

In the low temperature limit, we identified a universal behavior for various quantum 
gases. The thermodynamic curvature is a power-law function with respect to T. The 
exponent is related to the dimension of the system and the dispersion relation of the 
particles but it is independent of the details of the statistics. 
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